capture program drop get_mean
program define get_mean
	local varname =  e(depvar) 
	qui su `varname'
	estadd scalar av = round(r(mean),0.001), replace
end 
